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ABSTRACT 

We present the result of near-infrared (near-IR) [Fe II] line mapping of the 
supernova remnant IC 443 with the IRSF/SIRIUS, using the two narrow-band 
filters tuned for the [Fe II] 1.257 /im and [Fe II] 1.644 /iin lines. Covering a large 
area of 30' x 35', our observations reveal that [Fe II] filamentary structures exist 
all over the remnant, not only in an ionic shock shell, but also in a molecular 
shock shell and a central region inside the shells. With the two [Fe II] lines, we 
performed corrections for dust extinction to derive the intrinsic line intensities. 
We also obtained the intensities of thermal emission from the warm dust asso- 
ciated with IC 443, using the far- and mid-IR images taken with AKARI and 
Spitzer, respectively. As a result, we find that the [Fe II] line emission relative 
to the dust emission notably enhances in the inner central region. We discuss 
causes of the enhanced [Fe II] line emission, estimating the Fe + and dust masses. 

Subject headings: infrared: ISM -- ISM: individual objects (IC 443) -- ISM: 
supernova remnants 
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Introduction 



IC 443 is a Galactic supernova 



Pet 



c ct al. 



1988; 



Olbert et al. 



2001) 



remnant (SNR) with an age of 3000 —30000 years 



20061 ) and metal-rich X-ray plasma ( JTroja et al. 



. Existence of a p ulsar wind nebula (IGaensler et al. 



20081 ) indicates that the SNR is of a 



core-collapse origin. IC 443 shows a limb 



and radio wavelengths (IBykov et al. 



2008: 

' r 



brightened morphology at optical, infrared, 



Castelletti et al. 



centrally-peaked morphology in the X-ray ( jTroja et al. 



201l[ ). while it shows a 



20081 ). Hence IC 443 is recognized 



as a mixed-morphology -type SNR. The distance to IC 443 is estimated to be 1.5 kpc 



(I Welsh &: Sallmen 



20031 ) and we use this value throughout the present paper. 



One of the most interesting and impressive features of IC 443 is its heavy interaction 
with the interstellar medium (ISM) via ionic and molecular shocks, which makes a marked 
contrast between the near-infrared (near-IR) morphologies. The 2MASS survey revealed 
that IC 443 has a bright shell in the J and H bands in the northeast region, which is faint in 



the i f, band. From the shell, the [Fe II] 1.644 /im line emission was detected (JGraham et al. 



19871 ). 



Rho et al. 



( 120011 ) concluded that the dominant carriers of both J and H bands are 
[Fe II] line emission. In contrast, the southern ridge is bright in the K s band but faint in 
the J and H bands. The near-IR spectroscopic observation with AK ARI detected man y 



2011J ). 



ro-vibrational transition lines of H 2 from parts of the southern ridge (jShinn et al. 
There fore the H? emission a t 2.12 Mm is likely to dominate in the K s band in the southern 



ridge. 



Noriega-Crespo et al. 



(120091 ) showed that the [Fe II] 26 /mi line emission also exists 
in the spectra taken along the southern ridge with the Spitzer /IRS. Thus IC 443 is an 
excellent laboratory to investigate ionic and molecular shock interactions between a SNR 
and the ambient ISM. 

This paper focuses on the spatial distribution of the [Fe II] line emission in IC 443. 
The origin of the [Fe II] line emission in SNRs is still controversial. Some mechanism to 
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efficiently produce gas-phase Fe is required 



onto dust grains in the interstellar space ( iDraine 



because mor e than 99% of Fe is usually depleted 



19951 ). A conventional idea is that the 



gas-phase Fe is produced by sputtering Fe-bearing dust grains, but we cannot rule out a 
possibility of supernova ejecta origins. 



KooetaL 



(120071 ) suggested that the [Fe II] line 



emission in SNR Gil. 2-0. 3 is partly of an ejecta origin, because, considering its expansion 
velocity, an ejecta is likely to be excited recently by a reverse shock, where the [Fe II] line 
was detected. To examine the origin of the [Fe II] line emission, it is important to make 
a detailed comparison of the spatial distribution of the [Fe II] line with that of the dust 
emission associated with IC 443. 

2. Observations and data reduction 

We observed IC 443 with the near-IR camera SIRIUS (Simulta neous Infrared Imager 



for Unbiased Survey; 



Nagashima et al. 



1999 



Nagayama et al. 



20031 ) on the IRSF (Infrared 



Survey Facility) 1.4 m telescope. IRSF is located at the South African Astronomical 
Observatory (SAAO), and has been operated by Nagoya University and the SAAO. The 
SIRIUS camera has a field of view of 7'7 x I'.l with a pixel scale of 0'.'45. We observed 44 
contiguous fields toward IC 443, which, in total, correspond to the area of 30' x 35', using 
the two narrow-band filters tuned for the two [Fe II] lines at 1.257 /im and 1.644 /im. The 
details of the observations are summarized in Table [TJ 

We applied the standard data reduction procedure to the array images, including 
dark subtraction, flat-fielding, sky subtraction, and dithered-image-combining. As a 
result, we obtained 44 images for each of the 1.257 /im and 1.644 /im filters, which were 
combined into a large-area map, covering the main structure of IC 443. We performed 



photometric ca 



Skrutskie et al 



ibration, making comparison with the 2MASS point source catalog (PSC; 



20061 ). We used the J- and if -band magnitudes for the calibrations of the 
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1.257 /iin and 1.644 /im filter images, respectively, assuming that the magnitude of each 
star is the same between the corresponding band and filter images. For the photometry, we 
selected stars with fluxes higher than 12.5 mag and flux errors smaller than 0.05 mag from 
the 2MASS PSC, where we used sufficiently-isolated stars to avoid problems due to star 
crowding. The number of the stars used in the photometric calibration is typically 15 per 
image, and the uncertainties of the photo metric calibration coeffici ents are ~5 %. Finally, 



we subtract point sources with SExtracter (IBertin fe Arnouts 



19961 ) to avoid overestimating 



the [Fe II] line intensities due to the stellar flux. 

The far-IR images of IC 443 were derived from the AKARI all-sky-survey data. 
The AKARI all- sky-survey was performed between May 2006 and August 2007 in the 
cold mission phase with liquid helium cryogen, using the Far-Infrared Surveyor (FIS; 



Kawada et al. 



20071 ). The FIS has four far-IR photometric bands at central wavelengths of 



65 fjtm (N60), 90 fjtm (WIDE-S), 140 /im (WIDE-L), and 160 /an (N160), which have the 
effective band widths of 22, 38, 52, and 34 /mi, respectively. The data we re processed with 
the A KARI pipeline tool originally optimized for point source extraction (JYamamura et al. 



2009), and then with an additional pipeline to recover large-scale diffuse emission ( JDoi et al. 
20121 ). We also used the Spitzer /MIPS 24 /im image of IC 443, retrieved from the 
Spitzer archives flhttp: //archive. spitzer.caltech.edu). The observation of IC 443 was part 
of the LATT ERJ5 E program (PI: G. Rieke). The data were taken in a scan map mode 



(IRieke et al 



20041 ) and produced by the pipeline version S18.12.0. 



Results 



Figure [T] shows the images of IC 443 with the [Fe II] 1.257 //m and 1.644 /im filters, 
before the removal of the point sources. As can be seen in the figu re, we detect the [Fe II] 



filamentary structures from a wide area of IC 443. 



Rho et al. 



toll ) detected J- and H-btmd 
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diffuse emission with the surface brightnesses of 2.3 x 10 4 and 1.6 x 10 4 ergs s 4 cm 2 sr 1 , 
respectively, from a 48" x 18" area in the northe ast shell, centered at (R.A., Decl.) = 



(06 h 17 m 34!41, +22°52'55'.'2). At the same position. iGraham et all ( 119871 ) detected the [Fe II] 
1.644 /im line intensity of 1.1 x 10 -4 ergs s _1 cm -2 sr -1 , and our maps show the [Fe II] 
1.257 /im and 1.644 /im line intensities of 1.7 x 10~ 4 and 1.2 x 10 ergs s _1 cm -2 sr -1 , 
respectively. Thus in the northeast shell, indeed, the [Fe II] 1.257 /im and the 1.644 /im line 
emission are dominant in the J and H band, respectively. More importantly, the [Fe II] 
filaments are detected not only in the northeast shell, which is bright in the J and H 
bands, but also in other large area including t he sou thern ridge and a central region inside 



the shells. For the entire remnant, 



RhoetaL 



(1200 ll ) obtained the J- and if-band fluxes of 



1.2 x 10 9 and 6.4 x 10 10 ergs s 1 cm 2 , respectively, while we detected the [Fe II] 1.257 
/im and 1.644 /im line fluxes of 8.2 x 10~ 10 and 6.2 x 10~ 10 ergs s^ 1 cm -2 , respectively. 

To obtain the intrinsic [Fe II] line intensities, we performed corrections for foreground 
dust extinction, using the observed line ratio of the two [Fe II] lines at 1.257 /im and 1.644 
/im. Since they are due to elect ronic transitions from the sa me upper level, the intrinsic 



line ratio must be fixed at 1.36 ( Nussbaumer fe Storey 



19881 ). We estimated the foreground 



extinction by comparing the observed line r atio w ith the theoretical one, assuming the 



interstellar extinction law of 



Cardelli et al. 



(119891 ) . In calculating the line ratios, we first 



masked the pixels where [Fe II] 1.257 /im and 1.644 /im are below la, because the fluctuation 
of background can produce strong positive or negative spikes in a ratio map. Then we 
smoothed the resultant ratio map with a boxcar of 150 x 150 pixels. The extinction map 
thus derived is shown in the left panel of Fig. [5] on the color scale, with the contour s of the 
12 CO map integrated in the velocity range of —10 to km s _1 from lLee et al.l (120121 ) . From 
the figure, it can be seen that the foreground extinction is large in the west and southeast 



regions, and CO clouds are lying in similar regions where 



Troja et al 



(120061 ) showed the 



absorption of X-rays. Using the extinction map, we performed the extinction corrections 
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for both 1.257 /im and 1.644 jum images. The resultant [Fe II] 1.257 /im line intensity image 
is shown in the right panel of Fig. [2] on the gray scale. 

The far- and mid-IR images of IC 443, which are obtained by AKARI and Spitzer, 
respectively, include emission not only from the dust associated with IC 443, but also from 
foreground and background dust. Therefore discrimination of these two components is 
required. The dust associated with IC 443 is heated by shocks and therefore likely to have 
higher temperatures than foreground and background dust. We do not consider the dust 
component associated with IC 443, but not yet heated by shocks, if any, because it cannot 
be the origin of the observed Fe + through dust sputtering. Hence we fitted mid- to far-IR 
spectral energy distributions (SEDs) by a two-temperature modified blackbody model with 
the emissivity power- law index of unity. First, we fitted the model to the SED of the 
total region with a photometry radius of 19', centered at (R.A., Decl.) = (06 h 17 m 18!63, 
+22 o 37'07'.'5). As a result, we find that the model with the temperatures of T co id = 14.8 
K and T warm = 56.3 K reproduces the total SED very well. Then with the temperatures 
fixed at the above best-fit values, we fitted local SEDs for every spatial bin of 15" x 15", 
allowing only the amplitudes of the two components to vary. Here, considering possible 
severe contaminations of the [O I] 63 /iin line emission to the 65 /xm narrow-band (N60) 
intensity and [C II] 158 /im line emission to the 160 /zm narrow-band (N160) intensity, 
we excluded the N60 and N160 band intensities in the local SED fitting. The contours in 
the right panel of Fig. [2] show the resultant distribution of the intensity of the warm dust 
emission integrated over the wavelength range of 2 — 200 /im, which is derived from the 
above SED fitting. 

Here and hereafter we define regions A to D, as shown in the right panel of Fig. [2], all of 
which have the same area of 2.1 x 2.1 pc 2 (4'9 x 4(9). Regions A, centered at (R.A., Decl.) = 
(06 h 17 m 55!91, +22°46'24"8) and D at (R.A., Decl.) = (06 h 17 m 40:58, +22°23 / 55"0) represent 



the ionic shock and the molecular shock region, respectively, while regions B at (R.A., 
Decl.) = (06 h 17 m 13!62, +22°41 , 10'. , 0) and C at (R.A., Decl.) = (06 h 17 m 2r!75, +22 o 32'10 / . , 0) 
represent inner central regions where the [Fe II] line emission is detected significantly. 
As can be seen in the figure, both [Fe II] and warm dust emission are strong in region 
A, while only the warm dust emission is strong in region D. In regions B and C, the 
[Fe II] line emission is relati vely strong whereas the dust emission is faint. Near region C, 



Noriega-Crespo et al 



(120091 ) also detected the [Fe II] 26 /im emission in their Spitzer /IRS 
spectrum. More quantitatively, Fig. [3] shows the correlation plot between the intensities of 
the extinction-corrected [Fe II] 1.257 /zm line and the warm dust emission integrated over 
the wavelength range of 2 — 200 /im. For comparative purpose, we draw the line fitted to 
only the data points of region A, which gives the slope of 0.058 ± 0.006 and the correlation 
coefficient of 0.3. The figure clearly reveals that the [Fe II] line emission in regions B and C 
is notably strong as compared with the warm dust emission. 



4. Discussion 



We estimate the masses of Fe + and warm dust for regions A to D. In a two- level 
approximation, the [Fe II] 1.644 /xm line intensity is given by 

T4 -0.94 e -1.57/T 4(ArFe+/10 l 6) 



'A1.644 



1 + 4.2 x 10 4 T 4 a69 /n e ' 
where T 4 is the gas temperature in units of 10 4 K, and n e and Np e + are the electron 
densi ty and the Fe + column density in units of cm -3 and cm -2 , respectively (IBlietz et al. 



994T). Assuming T4 = 1.2 and n e = 500, which were derived from the mid-IR [Fe II] lines 



( JRho et al 



200ll), we estimate Np c + and thus the Fe mass, M Fe +, from the area of each 
region. As a result, M Fe + is 5, 1, 0.9, and 0.8 x 10~ 4 M for regions A to D. Note that, 
if we adopt the typical values derived from the near-IR [Fe II] lines in SNRs, T 4 = 0.5 



and n. 



1 x 10 4 fe. 



Graham et al. 



1987 : 



Koo et al. 



20071 ). M Fe + decreases by an order 
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of magnitude, whi 



e if we adopt t he gas density of only a few ti mes 10 cm 3 which was 



20081 ) and optical line observations (IFesen fc Kirshner 



19801 ). 



suggested by H I (ILee et al. 

M Fe + increases by an order of magnitude. Accordingly the above Fe + masses can have such 

systematic errors. 



On the other hand, by assuming the dust absorption coefficient given by 



Hildebrand 



( 119831 ). a grain radius of 0.1 /xm, and a specific dust mass density of 3 g cm 3 , the mass of 



the warm dust, Ma, is expressed as 



M d = 10" 



(- 



Mr: 



©) 



(2) 



\1Q 8 L Q J V40K 

where T warm and L VSLTia are the temperature and luminosity of the warm dust, respectively, 
derived from the above SED fitting. As a result, we obtain M<j of 50, 2, 0.9, and 50 
x 10~ 4 M Q for regions A to D. Hence the M Fe +/M d ratios are 0.1, 0.5, 1, and 0.02 for 
regions A to D. To first order, the absorption coefficient of a silicate grain is proportional 
to its volume, and thus the mass estimation is not affected by the grain size, although 



Andersen et al. 



the s hock processing of dust will change the grain size distribution (e.g. 

20111 ). If stochastic heating of very small grains is dominant, the dust masses could be 

overestimated. 

The above huge differences in M ¥e +/M d can be interpreted by the differences in the 
ratio of the destroyed to the non-destroyed dust mass from region to region. We estimate 
the latter ratios to be 30, 70, 80, and 8 % for regions A to D, assuming that all elements in 
dust are removed away by sputtering at the same rate, all gas-phase Fe is singly ionized, and 



Fe initially depleted to dust accounts for 22 % of M d ( JDraine et al.l 120071 ). Accordingly, the 
initial dust masses are 70, 7, 5 and 50 x 10 -4 M Q for regions A to D. The difference between 
regions A and D can be explained by the fact that more destructive shock is propagating in 
region A tha n in region D w here the deceleration of the shock due to the molecular cloud 



is indicated (IRho et al. 



20011 ) . Our result shows that the spatial distribution of the initial 
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dust grains is highly biased toward regions A and D, probably reflecting the distribution of 
the pre-existing ISM. 

In regions B and C, the dust seems to be almost completely destroyed, which might 
be reasonable considering that the shock front reached the interior regions earlier. The 
spu ttering destruction time of dust in the hot plasma of IC 443 is estimated to be ~1 x 10 5 



yr (IDraine fc Salpeter 
10 7 K (IPetre et al 



19791 ). by using the plasma density and temperature of 1 cm 3 and 
19881 ). respectively, and assuming the initial grain r adius of 0.1 nm. 



This time scale is, however, considerably longer than the age of IC 443. iJones et al 



(1994) 



showed that dust destruction associated with a J-shock is much more efficient, which is 
suitable for regions A and D, and could also be viable for regions B and C. On the other 



hand, the time for Fe to reach an ionizatio n equi 



n e is the electron density of X-ray plasma (Masai 



using n e of 1.7 cm 3 obtained for IC 443 by 



ibrium in hot plasma, ~10 12 n e x s, where 
19941). is est i mated to be ~2 x 10 4 yr, by 



Yamaguchi et al. 



(120091 ). Since this time scale 



is comparable to the age of IC 443, a significant fraction of gas-phase Fe atoms is likely to 
be more highly ionized than Fe + in regions B and C, whereby we may underestimate the 
Fe + masses (i.e., the destroyed dust masses) there. 

The Fe + mass in the central region of IC 443 is estimated to be 2 x 10~ 3 M from 
the total [Fe II] line flux in the observed area relatively free of the pre-existing ISM, where 
the intensity of the warm dust emission is lower than 1.5 x 10~ 3 ergs s _1 cm -2 sr _1 (i.e., 
the lowest contour level in the right panel of Fig. [2]). According to the model calculation, 
the 56 Fe mass produced from a 15 — 25 M & core-collapse supernova is 0.05 — 0.13 M 
( Thielemann et al.lll996l ). which is one to two orders of magnitude larger than the Fe + 



mass derived from our estimation. Thus our result suggests that Fe ejecta may not be so 
abundant as predicted theoretically for a core-collapse supernova. 
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Conclusion 



With the IRSF/SIRIUS, we observed the Galactic SNR IC 443, using the narrow-band 
filters for the [Fe II] 1.257 /zm and 1.644 /zm lines. Our observations detected the [Fe II] 
line emission from a wide area of IC 443, not only the shock-heated shell regions, but also 
interior regions, with prominent filamentary, shell-like structures. From a combination 
of the [Fe II] 1.257 /im and 1.644 /im lines which are attributed to electronic transitions 
from the same upper energy level, we estimated the foreground extinction and obtained 
extinction-corrected [Fe II] line maps. We also obtained the map of the shock-heated warm 
dust emission using AKARI and Spitzer IR images. As a result, we find that the [Fe II] 
line emission is notably strong relative to the warm dust emission in the inner central part 
of IC 443. We estimated the Fe + masses from the intensity of [Fe II] line and compared 
them with the masses of the shock-heated warm dust, which are derived from the SED 
fitting by using the AKARI and Spitzer IR images. The result implies that the mass ratios 
of Fe + to the dust are 0.1 (ionic shock shell), 0.02 (molecular shock shell), and 0.5 and 1 
(interior regions). The large difference in the ratio likely reflects that in the dust destruction 
efficiency; in particular, the Fe + /dust enhancement in the interior regions requires an 
efficient destruction mechanism for large grains, or contribution of theoretically-predicted 
Fe ejecta. 
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Table 1. Summary of the present IRSF observations 



-^center (/-^m) Transition AA e ff a (/xm) Exposure (s) NDI b Date 



1.257 
1.644 



a 4 D 7 / 2 ->■ a 6 D 9 / 2 



a 4 D 7/ o ->■ a 4 F ( 



y 7/2 



9/2 



0.028 
0.018 



60 
30 



10 



Nov 2012 



15 Feb, Nov 2012 



a Effective band width calculated from jS(X)d\ = T\AX e g, where S(A) is a filter re- 
sponse curve, and T A is a throughput at A cente r- 

b Number of dithered images. 
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Fig. 1.— [Fe II] 1.257 /iin (left) and 1.644 /im (right) line maps, smoothed with a Gaussian 
kernel of 2.3" in sigma. The color levels are given in units of ergs s _1 cm~ 2 sr _1 . 
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Fig. 2.— (Left) Foreground dust extinction Ay map derived from the observed [Fe II] line 
intensities, with th e contours of th e 12 CO map integrated in the velocity range of —10 to 



km s 1 taken from 



Leeetal 



( 120121 ). (Right) Spatial distribution of the extinction-corrected 
[Fe II] 1.257 pm line intensity, with the contours of the intensity of the warm dust emission 
integrated over the wavelength range of 2 — 200 //m, which is derived from the SED fitting. 
The gray scale ranges from 2.0 to 8.0 x 10~ 5 ergs s _1 cm~ 2 sr _1 , while the contour levels 
are 1.5, 2.5, and 3.5 x 10~ 3 ergs s _1 cm -2 sr _1 . 
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Fig. 3.— Correlation plot between the intensities of the extinction-corrected [Fe II] 1.257 
/im line and the warm dust emission integrated over the wavelength range of 2 — 200 /xm. 
Regions A to D are defined in Fig. |2} and the line fitted to only the data points of region A 
is overlaid. The data are spatially sampled every 1', which corresponds to 130 pixels of the 
camera array. 



